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ARTICLE INFO ABSTRACT 


Keywords: 
Pu-erh tea 


Moist-heat action and microbial action are essential to the formation of Pu-erh tea’s quality. To explore the 
impact of piling fermentation (PF) environment on Pu-erh tea quality, samples subjected to spontaneous PF 


Microbial action (SPPF) and sterile PF (STPF) from the same batch of raw material were processed, and their sensory quality and 


pay ig ‘ metabolic profile dynamics were determined. Sensory evaluation indicated that infusions of the SPPF samples 
Piling fermentation : : £ : 
Quality were mainly brownish-red, stale-mellow taste, and stale aroma, whereas infusions of the STPF samples were 


mainly orange-red, sweet-mellow taste, and woody aroma. Metabolomics analysis identified 465 key metabo- 
lites, including 33 microbial markers and 39 moist-heat activity markers. Enzyme and pathway analyses revealed 
that the production of microbial enzymes accelerates the degradation of flavonoids, flavonols, and lipids and also 
indicated an exclusive influence of microorganisms on the lipid metabolite pathway. This study elucidated the 


mechanism underlying the formation of Pu-erh tea’s stale flavour. 


1. Introduction 


Dark tea is a unique type of post-fermentation tea produced from the 
solid fermentation of tea leaves (Camellia sinensis) (Wang, Qiu, Gan, & 
Zhu, 2021; Zhu et al., 2020). On the basis of differences in processing 
and origin, dark tea is traditionally divided into Yunnan Pu-erh tea, 
Hunan Fuzhuan, Hubei Qingzhuan, and Guangxi Liubao tea (Lv, Zhang, 
Lin, & Liang, 2013). Pu-erh tea is produced in Yunnan Province, China, 
and can be classified into raw and ripe Pu-erh teas on the basis of the 
degree of fermentation (Cao et al., 2020; Wang et al., 2018). Piling 
fermentation (PF) is considered the process crucial to achieving the 
sensory attributes of Pu-erh tea (brownish-red tea infusion, stale and 
fungal aroma, and stale-mellow taste) and is called the “quality flavour 
process” (Cao et al., 2020; Wang et al., 2018; ISO, 2023; Zhang, Li, Ma, 
& Tu, 2011). As mentioned in our previous study, three theories describe 
the PF of dark tea: enzyme regeneration theory, microbial theory, and 
moist-heat theory. We considered it meaningful to continue to investi- 
gate how microbial action and hygrothermal action contribute to the 
change in tea quality during fermentation, building on the foundation 


laid by other researchers. 

Scholars have explored the transformation of many types of chemical 
components of dark tea. In reviews, Lv et al. and Zhang, Zhang, Zhou, 
Ling, & Wan (2013) suggested that the special sensory characteristics of 
Pu-erh tea arise from the various chemical changes and transformations 
of the chemical constituents of the sun-dried green tea (SDGT) leaves 
that occur during the PF process (Lv et al., 2013; Zhang et al., 2013). 
Many functional components of Pu-erh tea have been isolated and 
identified. Pu-erh tea is becoming increasingly popular, and researchers 
have developed a keen interest in the formation and changing quality 
patterns of Pu-erh tea. Researchers have reported that microorganisms 
are crucial to the formation of high-quality Pu-erh tea (Wang et al., 
2021). With the gradually increasing popularity of Pu-erh tea globally, 
research on Pu-erh tea is rapidly increasing, and a particular boom 
occurred between 2015 and 2022. In 2016-2017, Lin et al. and Wan 
et al. investigated the chemical profile, aroma profile, and antioxidant 
activity of finished dark teas obtained through various processes (Cao 
et al., 2016; Lv, Zhang, Shi, & Lin, 2017). Cheng et al. (2021) employed 
nontargeted metabolomics and electronic tongue detection to determine 
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the chemical changes and taste characteristics of five typical dark teas, 
providing a clear understanding of the flavour characteristics of Pu-erh 
tea and the differences between flavour and other tea qualities. Through 
computation of odour activity values and multivariate analysis, Pang 
et al. (2019) compared the potent active odour components in raw and 
ripe Pu-erh tea. Cao et al. (2020) used a metabolomics approach based 
on ultra-high-performance liquid chromatography (UHPLC) triple 
quadrupole mass spectrometry (MS) to systematically identify bio- 
markers of and various metabolites in raw and ripe Pu-erh tea. These 
studies have provided an accurate understanding of the effect of PF on 
the flavour of Pu-erh tea. In 2020-2021, Ning et al. thoroughly inves- 
tigated the aroma profile and flavour features of and the corresponding 
volatile and nonvolatile compounds in Pu-erh tea stored for various 
lengths of time (Xu et al., 2019; Xu et al., 2021), revealing how me- 
tabolites are affected by storage time. Shi et al. (2021) produced Pu-erh, 
Liubao, Qingzhuan, and Fuzhuan teas by using a single raw tea material 
and applied widely targeted metabolomics and gas chromatography—MS 
to investigate the impacts of various microbial fermentation methods on 
the teas’ chemical and volatile profiles. Deng et al. analysed the evolu- 
tion of the active flavour components of Pu-erh tea during artificial 
fermentation and established the tea’s flavour profile and metabolite 
change pattern (Deng et al., 2021). Researchers have mostly focused on 
the effects of processing, storage environment, tea type, and the tea 
leave themselves on the quality of Pu-erh tea. Investigations into the 
effects of microbial action and moist heat in the fermentation environ- 
ment on the quality of tea are rare. 

In our previous study, we employed metabolomics to compare raw 
SDGT samples with finished samples obtained through spontaneous PF 
(SPPF) and sterile PF (STPF) (Li et al., 2022). Through metabolomic 
analysis, we could screen for metabolites such as quercetin, dulcimer, 
and kaempferol which are potential markers of microbial action, and for 
ascorbic acid, $-glucuronide, and ellagic acid, which are potential 
markers of hygrothermal action. We discovered that kaempferol, quer- 
cetin, and ellagic acid were the main active substances formed under 
microbial action (Li et al., 2022). Both microbial and moist-heat actions 
were discovered to have strong effects on the quality of Pu-erh tea. 
However, how these marker metabolites change during processing and 
what factors regulate the synthesis and degradation of metabolites 
during fermentation; currently, remains unclear. Therefore, in this 
study, different assay techniques (non-targeted metabolomics, quanti- 
tative assay techniques, different enzyme activity assays, electronic 
tongue, etc.) were used for the first time to detect and analyze the 
non-volatile components of SPPF and STPF processes of Pu-erh tea. 
Combining multivariate statistical analysis to deeply investigate the 
influence of PF process on the quality formation mechanism of Pu-erh 
tea. The results of the present study provide a better understanding of 
the influence of SPPF and SPPF on the metabolic profile of Pu-erh tea 
and thus of the mechanism by which the quality characteristics of Pu-erh 
tea develop. 


2. Materials and methods 
2.1. Pu-erh sample preparation 


Samples were prepared at Menghai Tea Factory. And the preparation 
process was consistent with our previous study (Li et al., 2022). In brief, 
in the pre-PF stage, we mixed 18 tons of SDGT raw materials with water 
at a 3:1 mixing ratio. The tea pile was turned over about every seven 
days and samples were collected and stored in —80 refrigerator. A total 
of five pile turnings were performed. The following SDGT samples were 
collected: samples of SDGT processed through SPPF for 0 (SPPF-1), 7 
(SPPF-2), 14 (SPPF-3), 22 (SPPF-4), 32 (SPPF-5), and 38 (SPPF-6) days 
and samples of SDGT processed through STPF for 0 (STPF-A), 7 
(STPE-B), 14 (STPF-C), 22 (STPF-D), 32 (STPF-E), and 38 (STPE-F) days 
(Graphical abstract). 
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2.2. Chemicals and reagents 


Chemical reagents were purchased from Yuanye Bio-Technology 
(Shanghai, China), Sigma-Aldrich Fluka (St. Louis, MO, USA), Abbki- 
nev and other companies. Detailed information is provided in the sup- 
plementary materials S1. 


2.3. Sensory evaluation 


The tea samples were assessed by six experts (Li et al., 2022). To 
ensure unbiased evaluation, the tea samples were assigned random 
numbers and prepared by steeping 3 g of tea leaves in 150 mL of puri- 
fied, boiled water in a white porcelain cup for 5 min. Then the infusion 
was transferred to a white porcelain bowl for evaluated by GB/T, 2017 
standards (GB/T, 2017). 


2.4. Nontargeted metabolomic analysis 


In this study, we used our previous research methods 
UHPLC-Orbitrap—MS method. Briefly, first, 0.4000 + 0.0010 g of each 
sample was weighed in a 10 ml centrifuge tube; then, 8 ml of methanol/ 
water (7/3 v/v) was added, containing the internal standard at a con- 
centration of 5 mg/L. Ultrasonication (3500 Hz, 30 min/per session) was 
performed at room temperature at 4 h intervals, twice in total. Subse- 
quently, 5 ml of the supernatant was taken and centrifuged; 100 p11 was 
then aspirated for a 40-fold dilution. Then, it was placed in a liquid 
phase vial until LC-MS analysis was performed. And detailed informa- 
tion is provided in the supplementary materials $1 and S2. (Li et al., 
2022; Table S1). 


2.5. Qualitative and quantitative detection of quality components 


Different catechins and different amino acids in the samples were 
assayed by HPLC and high-speed amino acid analyzer respectively. 
These assays are consistent with our previous assays (Li et al., 2022; 
Ning et al., 2016). Detailed materials are provided in Supplementary 
Material S1. 


2.6. Extraction of metabolomics data 


To initiate the analysis of the LC-MS data set, the data were con- 
verted from raw format to mz/ML format and analysis base file format 
using MS-Convert and the Analysis Base File Converter. The data were 
processed by noise filtering, peak alignment, and peak filling in MS-DIAL 
version 3.82; secondly, the data were identified based on MS-FINDER 
version 3.04. In MS-DIAL version 3.82, noise filtering, peak alignment, 
and peak filling were performed on the data; secondly, the data were 
characterized based on MS-FINDER version 3.04, which mainly included 
m/z values, retention times, and relative MS intensities. Times, and 
relative MS intensities (Tsugawa et al., 2015). These software programs 
enable the incorporation of various databases to identify unknown 
compounds, making them versatile tools for metabolite identification. 
To ensure the data were of high quality, ions with a filling rate of less 
than 33% or a signal/noise <10 were deleted, and metabolites with 
relative standard deviations of less than 30% were analysed further. The 
identification of compounds was carried out by searching for metabo- 
lites in two databases: PubChem and Tea Metabolome (Li et al., 2022). 


2.7. Key enzyme activity assays 


The study analysed the activities of various enzyme, specifically 
anthocyanin reductase (ANR), polyphenol oxidase (PPO), flavonol syn- 
thase (FLS), pectase, peroxidase (POD), pectase, anthocyanidin synthase 
(ANS), phosphatidylcholine diglyceride choline phosphotransferase 
(PDCT), flavanone-3-hydroxylase (F3H), and Icithin acyltransferase 
(LPCAT). In brief, the enzyme activity assay analysis method involves 
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reagent preparation and an assay procedure. Reagents, including a wash 
solution, sample diluent, HRP-conjugate reagent, chromogen solution A, 
chromogen solution B, and stop solution, must be prepared before 
beginning the assay. The standard wells and testing sample wells are 
prepared by adding standards and samples in duplicate to a Micro Elisa 
Strip plate. The standard well is set and the testing sample well is added 
with a mixture of 10 yl testing sample and 40 pl sample diluent, while 
the blank well is left empty. The assay procedure includes adding 50 ul 
of standard to the standard well and 100 pl of HRP-conjugate reagent to 
each well, which is then covered with an adhesive strip and incubated 
for 60 min at 37 °C. The wells are then aspirated and washed with wash 
solution, and the chromogen solutions A and B are added to each well. 
After incubating for 15 min at 37 °C, the stop solution is added to each 
well, and the optical density (O.D.) is read at 450 nm using a microtiter 
plate reader within 15 min. The standard curve is then used to determine 
the amount in an unknown sample. The supplementary material pro- 
vides detailed information on the methodologies used to assay these 
enzyme activities (supplementary materials S1). 


2.8. Statistical analysis 


SIMCA-P 14.1 software was employed to conduct principal compo- 
nent analysis (PCA), and hierarchical cluster analysis (HCA) on the 
metabolomics data. Model quality was assessed using R? and Q”. To 
ensure the analysis was performed on only the most significant metab- 
olites, a rigorous statistical approach was employed. Specifically, me- 
tabolites were selected for further analysis (p < 0.05, VIP > 2, 1.50 < 
fold change and fold change <0.50), indicating a significant difference. 
Conversely, metabolites with p values greater than 0.05 were excluded 
from the analysis to prevent false positives and ensure the integrity of 
the results. To gain insight into the complex relationships between 
metabolites, Spearman analysis was used. And SPSS 25.0, TBtools, and 
Origin (version 2022) were performed. 


3. Results and discussion 
3.1. Analysis of sensory quality 


The SPPF- and SPPF-processed samples from the leaves of SDGT were 
evaluated in terms of the colour of their infusion colour, their taste, and 
their aroma. The results showed that the two fermentation environments 
resulted in Pu-erh tea of significantly different quality. As shown in 
Table S2, we found that the highest quality was achieved through SPPF 
(under microbial and moist-heat action) for up to 38 days; the infusion 
colour was brownish-red, and the aroma was dominated by notes of 
staleness, fungal staleness, caramel, and purity; the taste was dominated 
by staleness, mellowness, and purity. Regarding STPF, the tea leaves are 
mainly affected by moist-heat action during the STPF process. The 
quality characteristics of the STPF tea leaves were an orange-red and 
dull-red infusion colour, a woody aroma, a sweet aroma, and a sweet and 
mellow taste: lose the unique quality characteristics of Pu-erh tea. The 
aroma is stale aroma and fungus stale, the taste is stale and mellow, pure 
and normal, and the soup color is brownish-red are the unique quality 
characteristics of Pu-erh tea. We believe that the microbial action was 
the cause of this quality alteration. How changes in the processing 
environment during processing regulate the synthesis of metabolites was 
still unclear, however. Therefore, for further investigation of the mi- 
crobial and moist-heat effects on the quality of Pu-erh tea, nontargeted 
metabolomics was performed to examine the metabolites generated 
during processing and thus investigate the causes of Pu-erh tea quality 
formation. 


3.2. Analysis of key quality components in Pu-erh tea 


In Pu-erh tea processing, the oxidation and degradation of amino 
acids and catechins have a significant effect on the aroma and taste of tea 
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leaves (Wang et al., 2022). The analysis of key quality components 
revealed that large amounts of amino acids and catechin monomers 
were degraded during the leaves’ fermentation (Figures $1, S2 and S3). 
Under SPPF, the number of ester-type catechins (EGCG and ECG) 
decreased dramatically (cliff-like decline) from the first stage of 
fermentation, whereas the number of simple catechins (EGC and EC) 
increased in the second stage of fermentation, after which it rapidly 
declined (Fig. S1). Combining these findings with the STPF analysis, we 
concluded that in the first stage of fermentation, simple catechins are 
produced through the rapid degradation of ester-type catechins under 
the action of mainly microorganisms (Fig. $1), which results in more 
production than degradation of simple catechins, leading to an increase 
in simple catechin content. We also discovered that the tea’s GA content 
increased rapidly under microbial action, whereas it increased linearly 
under the action of moist heat, and eventually, the GA tended to 
degrade. The different amino acids examined are detailed in Figs. S2 and 
S3. Under natural fermentation, serine (Ser), aspartic acid (Asp), argi- 
nine (Arg), tyrosine (Tyr), phenylalanine (Phe), threonine (Thr), Glu- 
tamic acid (Glu), proline (Pro), Alanine (Ala), B-Ala, and theanine were 
degraded mainly under the microbial influence, whereas ornithine 
(Orn), glycine (Gly), and Cysthi were produced. Additionally, we found 
that Ile, Hylys, Val, and theanine were largely degraded by moist heat, 
whereas Orn was produced. These results yielded new insight into the 
mechanism of degradation of Pu-erh tea metabolites and a new under- 
standing of the formation of tea quality. 


3.3. Nontargeted metabolomic analysis 


3.3.1. Comprehensive metabolite analysis 

We extracted 22,814 (negative ion modes) and 10,256 (positive ion 
modes) metabolite features, respectively, from the 13 groups of samples. 
After data preprocessing, 6171 (negative ion modes) and 6725 (positive 
ion modes) metabolites identified through the positive and negative ion 
modes were used for analysis. Fig. 1 shows the results of performing PCA 
to compare the metabolites present in the samples. The PCA analysis for 
the positive ion mode (Fig. 1a) showed that the first two principal 
components, PC1 and PC2, accounted for 35.5% and 15.2% of the total 
variance, respectively (R?X = 0.653). The samples were clearly distin- 
guished from each other on the PCA score plot (Fig. 1a). In addition, the 
positioning of the quality control samples at the origin of the coordinates 
provided evidence of the accuracy of our detection data. An HCA 
dendrogram was produced based on the PCA findings to aid in the 
visualization of the variations between the samples. The tea samples 
were separated into three groups at a Euclidean distance of around 
5200-6000:raw (SDGT, SPPF-1, STPF-A), SPPF-processed (SPPF-2 to 
SPPF-6), and STPF-processed (STPF-B to STPF-E) samples (Fig. 1d). And 
the HCA dendrogram clearly distinguished the SDGT samples from the 
other three sample types, indicating that the chemical composition of 
the SDGT samples differed from that of the other samples. Regarding the 
negative ion mode (Fig. 1e—-h), on the basis of the PCA (PC1 = 0.29%, 
PC2 = 0.15%, R°X = 0.835) and HCA analysis, we reached the same 
conclusion—that the raw material could easily be distinguished from the 
SPPF and STPF samples. 

PLS-DA provides a more robust supervised classification method 
than PCA, allowing for better noise reduction and more valuable infor- 
mation extraction. In order to identify the metabolites impacted by SPPF 
and STPF, we performed PLS-DA and generated S-plots (Fig. 1b, c, f, and 
g). Based on the PLS-DA findings, an S-plot was constructed, and me- 
tabolites were extracted (Fig. 1d-f, VIP > 2). Each PLS-DA involved 
seven-fold cross validation with 50 permutations (Fig. 1c and g). The 
intercepts of R? and Q? for the positive and negative models depicted in 
Fig. 1¢ [R? = (0.00, 0.30), Q? = (0.00, —0.40)] and Fig. 1g [R* = (0.00, 
0.80), Q? = (0.00, —0.03)] indicated that the original model was 
feasible. The results of the cross-validation indicated that the metabo- 
lites were significantly impacted by SPPF and STPF. Statistical analysis 
identified a total of 206 (positive ion modes) and 259 (negative ion 
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Fig. 1. Metabolic profile analysis. (a and e) Score plots of principle component analysis. (b and f) Partial least squares—discriminant analysis. (c and g) S-plot analysis. 


(d and h) Score plots of hierarchical clustering analysis. 


Notes: quality control (QC), sun-dried green tea (SDGT), spontaneous piling fermentation (SPPF), and sterile piling fermentation (STPF). (a—d) Negative ion mode and 
(e-h) positive ion mode. (For interpretation of the references to colour in this figure legend, the reader is referred to the Web version of this article.) 


modes) metabolites (Fig. 2 and Fig. S4), respectively, with VIP values 
greater than 2 and p-values less than 0.05 (Fig. 3 and Fig. S5). As pre- 
sented in Table 1, the structures of the 465 metabolites were tentatively 
determined. 


3.3.2. Metabolite profile analysis 

Heat maps were generated to visualise the metabolite patterns in 
different ion modes. Fig. 2 and Fig. S4 shows the pattern of metabolite 
changes in the two modes, with these changes divided into five main 
categories. As illustrated in Fig. 2, for the negative ion mode, the 
following was discovered: Cluster I comprised metabolites mainly 
influenced by microbial action; cluster II comprised metabolites mainly 
degraded by the action of moist heat; cluster III comprised metabolites 


degraded by the combined effect of microbial action and moist heat; 
cluster IV comprised metabolites mainly degraded through microbial 
action, with the effect of moist heat being relatively small; and cluster V 
comprised metabolites synthesised through the effect of humid heat. 
However, for the positive ion mode (Fig. S4), the following was deter- 
mined: Cluster I comprised metabolites degraded mainly through mi- 
crobial and moist-heat actions; cluster II comprised metabolites 
degraded mainly through microbial action, with moist heat having little 
effect; cluster IV comprised metabolites that were not strongly influ- 
enced by their environment; cluster III comprised metabolites produced 
under the action of microorganisms, with some relatively unaffected by 
the actions of heat and humidity; and cluster V comprised metabolites 
produced under the moist-heat action, with possible degradation 
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Fig. 2. Heat map resulting from nontargeted metabolomic analysis: negative ion modes. 


through microbial action. To the best of our knowledge, we are the first 
to investigate the effects of SPPF and STPF on the quality of Pu-erh tea 
and to comprehensively analyze the effects of PF processes on the 
metabolic profile. 


3.3.3. Differential metabolite screening 

To further screen for key metabolites affected by the processing 
environment, we used volcano plots [| log2 (fold change) | > 1, —log10 
(P)] and Venn diagrams for in-depth analysis of the 465 metabolites 
(Fig. 3 and Fig. S5). As illustrated in Figs. S5a and S5b we could conclude 
that 36 of the metabolites were affected by the combined action of mi- 
crobes and moist heat; 10 (1 + 9) were downregulated, and 26 (11 + 15) 


were upregulated. As revealed in Table 1, the 36 metabolites affected by 
the action of moist heat and microbes included six flavonoids (kaemp- 
ferol, quercetin, puerin ii, puerin iii, puerin iv, and puerin i), six orga- 
nooxygen compounds [sphinganine, hexadecasphinganine (1+), 2- 
amino-1,3,4,5-icosanetetrol, N-oxide dodecyldimethylamine, p-tartaric 
acid, and shikimic acid], and three pyrimidine nucleosides. We found 
that most of these metabolites were synthesised through the action of 
moist heat and microbes. Further, 33 metabolites were affected by mi- 
crobial action, including eight (4 + 4) that were downregulated and 25 
(19 + 6) that were upregulated (Figs. S5a and S5b). As shown in Table 1, 
among the 33 metabolites influenced by microorganisms, there were 
nine glycerophospholipids [LysoPE(0:0/18:2(9Z,12Z)), LysoPC(18:2 
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Fig. 3. Differential metabolite screening (volcano plots). Note: sun-dried green tea (SDGT), spontaneous piling fermentation (SPPF), and sterile piling fermentation 
(STPF). (For interpretation of the references to colour in this figure legend, the reader is referred to the Web version of this article.) 


(9Z,12Z)), and LysoPE(18:2(9Z,12Z)/0:0)) etc.], five organo-oxygen 
compounds (choline, galactitol, arabinonic acid, p-glucose, 3-hydroxyir- 
isquinone, and 5-hydroxyirisquinone), four carboxylic acids (cis-diho- 
moaconitic acid, N2-succinyl-L-ornithine, 4-glutamylamino butanoate, 
and rosmarinic acid), and some derivatives. We discovered that these 
metabolites were mainly synthesised under the action of microorgan- 
isms; therefore, we concluded that this was the cause of the final quality 
formation of Pu-erh tea. Finally, 39 metabolites were affected only by 
moist-heat action, 17 (8 + 9) which were downregulated and 22 (13 + 
9) which were upregulated (Figs. S5a and S5b). As shown in Table 1, 
these 39 metabolites included 6 flavonoids [(—)-catechin 3-O-gallate, 
dihydroisorhamnetin, apigenin 7-xyloside, apigenin 7-3’-acetyl-6’-E-p- 
coumaroylglucoside, epicatechin quinone, and dihydrokaempferol], 10 
organo-oxygen compounds (phytosphingosine, phytosphingosine, 
quinic acid, serotinose, serotinose, quinic acid, linalool 3,6-oxide pri- 
meveroside, N-acetylneuraminic acid, phenylethyl primeveroside, and 
arbutin), 3 carboxylic acids (UNPD127787, acetylenedicarboxylate, and 
2-piperidinecarboxylic acid), and some derivatives. We found that most 
of these metabolites were degraded under the action of humid heat and 
thus concluded that the degradation of metabolites was promoted under 
the action of humid heat. In summary, we screened a total of 108 me- 
tabolites that were significantly affected by microbial and moist-heat 
effects and determined that these metabolites are vital to the quality 
formation of Pu-erh tea (Table 1). 


3.4. Effects of SPPF and STPF on the metabolic pathways of Pu-erh tea 


3.4.1. Kyoto encyclopedia of genes and genomes metabolic pathway 
analysis of various metabolites 

On the basis of the 108 key metabolites obtained from the screening, 
we performed KEGG pathway analysis to screen the key pathways 
affected by microbial and moist-heat actions (Fig. S6). KEGG metabolic 
pathway analysis is useful for understanding the functional information 
of secondary metabolites in metabolic pathways (Kumari & Parida, 
2018). From the KEGG database, 21 relevant pathways covering 55 
metabolites were identified (Fig. S6). We found that the two most 
relevant metabolic pathways were the sphingolipid and glyceropholipid 
metabolic pathway (mainly including glyceropholipid, glyceropholipid, 
and glyceropholipid; 18 metabolites in total) and the flavone, flavanols, 
and flavonoid biosynthesis pathway (mainly including quercitrin and 
quercetin, and 10 metabolites in total; Fig. 4). To clearly visualise the 
metabolite changes in the aforementioned two metabolic pathways, we 
created a metabolic pathway by querying the aforementioned 
KEGG-identified metabolic pathways. By consulting previous studies, 
we discovered that the degradation and synthesis of flavonoids play vital 
roles in the quality (aroma and taste) of Pu-erh tea (Li et al., 2022). 
However, the phospholipid pathway is currently little studied, and to 
our best knowledge, there is only one relevant study in tea. 

To better understand the mechanism through which metabolites are 
regulated during fermentation, the enzyme activities and metabolite 
contents corresponding to the relevant pathways were measured and 
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Table 1 
Metabolite identification. 
Factors Down/Up m/z Formula Adduct ions Metabolite name Compounds ID QLs 
classification 
Microbial Down 137.0244 C7H603 M-H]- 4-Hydroxybenzoic acid Benzene and 1234 b 
action (33) regulation substituted derivatives 
(8) 361.0893 C18H1608 M+H] Rosmarinic acid Cinnamic acids and 5617 c 
derivatives 
478.2928 C23H44NO7P M+H] LysoPE(18:2(9Z,12Z)/0:0) Glycerophospholipids 7622 c 
564.3303 C35H43N502 M-H]- UNPD194693 Harmala alkaloids 17240 c 
721.3647 C42H50N407 M-H]- UNPD169811 Indoles and derivatives 20532 c 
165.0402 C5H1006 M-H]- Arabinonic acid Organooxygen 2291 ( 
compounds 
460.2691 C26H37NO6 M+H] Militarinone A Prenol lipids 7346 c 
415.2114 C24H3006 M+H] Eplerenone Steroids and steroid 6542 ¢c 
derivatives 
Up 171.0288 | C7H605 M+H] Gallic acid Benzene and 1682 a 
regulation substituted derivatives 
(25) 203.0527 C8H1006 M+H] cis-dihomoaconitic acid Carboxylic acids and 2366 c 
derivatives 
233.1132 C9H16N205 M+H] N2-Succinyl-l-ornithine Carboxylic acids and 2954 c 
derivatives 
233.1132 C9H16N205 M+H] 4-(Glutamylamino) butanoate Carboxylic acids and 2953 c 
derivatives 
112.0505 C4HS5N30 M+H] Cytosine Diazines 192 b 
277.0893 C12H18N203 M+K] Secobarbital Diazines 3810 c 
365.1054 C12H22011 M+NH4]+ Sucrose Disaccharides 5669 b 
295.2275 C18H3203 M-H]- Alpha-dimorphecolic acid Fatty Acyls 7360 ( 
355.1049 C16H22010 M-H20-H]- Geniposidic acid Flavonoids 9657 c 
562.3145 C25H43N0O10 M-+FA-H]- LysoPC(18:3(6Z,9Z,12Z)) Glycerophospholipids 17163 (a 
564.33 C26H50NO7P M-+FA-H]- LysoPC(18:2(9Z,12Z)) Glycerophospholipids 17239 c 
476.2771 C23H42NO7P M+H] LysoPE(0:0/18:3(6Z,9Z,12Z)) Glycerophospholipids 7600 c 
518.3242 C26H48NO7P M+H] LysoPC(18:3(6Z,9Z,12Z)) Glycerophospholipids 8107 c 
522.3555  C26H52NO7P M+H] LysoPC(18:1(9Z)) Glycerophospholipids 8145 c 
520.3397 C26H50NO7P M+H] LysoPC(18:2(9Z,12Z)) Glycerophospholipids 8126 c 
476.2772 C23H42NO7P M+H] LysoPE(0:0/18:3(9Z,12Z,15Z)) Glycerophospholipids 7601 c 
478.2929 C23H44NO7P M+H] LysoPE(0:0/18:2(9Z, 12Z)) Glycerophospholipids 7623 c 
206.0457. CSH13NO4S M+Na]+ Choline sulfate Organic sulfuric acids 2408 c 
104.107 CSH14NO M]+ Choline Organonitrogen 53 c 
compounds 
217.0482 C6H1406 M+Cl]- Galactitol Organooxygen 4522 c 
compounds 
181.0721 C6H1206 M+H] d-Glucose Organooxygen 1935 c 
compounds 
391.2843 C24H3804 M+H] 3-Hydroxyirisquinone;5- Organooxygen 6104 c 
Hydroxyirisquinone compounds 
266.0746 C9H13N305 M+Na]+ Cytidine Pyrimidine nucleosides 3585 c 
474.2622 | C30H37NO4 M-H]- ulipristal acetate Steroids and steroid 14298 c¢ 
derivatives 
609.2706 C35H36N406 M+H] epoxypheophorbide a Tetrapyrroles and 8968 c 
derivatives 
Microbial Down 149.0234 C8H403 M+H] phthalic anhydride Benzofurans 1170 c 
action and regulation 174.956 C3H603S2 M+Na-2H]- Bruguiesulfurol Dithiolanes 2705 c 
Moist-heat (10) 256.2635 C16H33NO M+H] Palmitic amide Fatty Acyls 3375 c 
action (36) 341.3049 C21H4003 M+H] Polyoxyethylene (600) monoricinoleate Fatty Acyls 5243 c 
282.279 C18H320 M+NH4]+ (9Z,12Z,15Z)-octadecatrien-1-ol Fatty Acyls 3932 c 
302.3053 C18H39NO2 M+H] Sphinganine Organonitrogen 4373 c 
compounds 
274.2739 C16H35NO2 M+H] hexadecasphinganine(1+) Organonitrogen 3754 C 
compounds 
362.3262 C20H43NO04 M+H] 2-amino-1,3,4,5-icosanetetrol Organonitrogen 5639 c 
compounds 
230.2478 C14H31NO M+H] dodecyldimethylamine N-oxide Organonitrogen 2912 c 
compounds 
437.1932 C24H3006 M+Na]+ Eplerenone Steroids and steroid 6906 c 
derivatives 
Up 400.1397 C21H23NO7 M-H]- Clivacetine Amaryllidaceae 11254 ¢ 
regulation alkaloids 
(26) 400.1399 C21H23NO7 M-H]- Clivacetine Amaryllidaceae 11257 [e 
alkaloids 
570.1606 C28H27NO12 M+H] (+)-Lactonamycin Anthracenes 8592 c 
570.1607 C28H27NO12 M+H] (+)-Lactonamycin Anthracenes 8601 c 
207.0509 C8H8N403 M-H]- Aminofurantoin Azolidines 4114 c 
468.127 C21H27NO9S M-H]- 3-(3,4-dimethoxyphenyl)-N-[2-(3,4- Benzene and 14060 c 


dimethoxypheny])ethyl]-3-(sulfooxy) substituted derivatives 
propanimidic acid 


(continued on next page) 
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Factors Down/Up m/z Formula Adduct ions Metabolite name Compounds ID QLs 
classification 
136.0618 C4H604 M+NH4]+ Succinic acid Carboxylic acids and 743 b 
derivatives 
554.1656 C28H27NO11 M+H]+ Amidepsine B Depsides and 8448 c 
depsidones 
175.0246 C6H806 M-H]- Ascorbic acid Dihydrofurans 2747 b 
285.0402 C15H1006 M-H]- Kaempferol Flavonoids 6940 a 
301.035 C15H1007 M-H]- Quercetin Flavonoids 7557 a 
402.1548 C21H23NO7 M+H] Puerin II Flavonoids 6304 b 
402.1549 C21H23NO7 M+H] Puerin III Flavonoids 6306 b 
402.1549 C21H23NO7 M+H] Puerin IV Flavonoids 6307 b 
402.1548 C21H23NO7 M+H] Puerin I Flavonoids 6300 b 
152.0567 CSHS5N50 M+H] Guanine Imidazopyrimidines 1228 b 
136.0618 = CSH5N5 M+H] Adenine Imidazopyrimidines 742. b 
149.009 C4H606 M-H]- d-Tartaric acid Organooxygen 1663 c 
compounds 
175.0577. C7H1005 M+H] Shikimic acid Organooxygen 1819 b 
compounds 
293.2473 C18H280 M+CH30H+H]+ __ Dinortrixagone Prenol lipids 4206 c 
268.104 C10H13N504 M+H] Adenosine Purine nucleosides 3619 c 
290.0859 C10H13N504 M+Na]+ Adenosine Pyrimidine nucleosides 4131 c 
244.0926 C9H13N305 M+H] Cytidine Pyrimidine nucleosides 3132 c 
400.1399 C21H23NO7 M-H]- Margrapine A Quinolines and 11255 oc 
derivatives 
555.2841 C29H4408 M-Cl]- Digoxigenin monodigitoxoside Steroids and steroid 17014 ~¢ 
derivatives 
331.0669 C13H16010 M-H]- beta-Glucogallin Tannins 8772 b 
Moist-heat Down 139.039 C7H603 M+H] Salicylic acid Benzene and 809 b 
action (39) regulation substituted derivatives 
(17) 112.9857 C4H204 M-H]- acetylenedicarboxylate Carboxylic acids and 229 c 
derivatives 
114.0915 C6H11NO M+H] 2-Piperidinecarboxylic acid Carboxylic acids and 281 b 
derivatives 
153.0182 C4H802S2 M+K] Oxidized dithiothreitol Dithianes 1261 c 
317.0542 C15H1008 M-H]- Dihydroisorhamnetin Flavonoids 8302 c 
401.0873 C20H1809 M-H]- Apigenin 7-xyloside Flavonoids 11301 ¢ 
619.1458 C32H28013 M-H]- Apigenin 7-(3’-acetyl-6"-E-p- Flavonoids 18846 c¢ 
coumaroylglucoside) 
289.0706 C15H1206 M+H] Epicatechin quinone Flavonoids 4111 
289.0706 C15H1206 M+H] Dihydrokaempferol Flavonoids 4105 
195.0507. C7H8N403 M-H]- 1,7-Dimethyluric acid Hydroxy acids and 3640 c 
derivatives 
191.056 C7H1206 M-H]- Quinic acid Organooxygen 3475 b 
compounds 
311.0981 C11H20010 M-H]- Serotinose Organooxygen 7965 c 
compounds 
348.0702 C11H19NO9 M+K] N-Acetylneuraminic acid Organooxygen 5352 c 
compounds 
439.1573, C19H28010 M+Na]+ Phenylethy! primeveroside Organooxygen 6961 b 
compounds 
311.0523 C12H1607 M+K] Arbutin Organooxygen 4577 c 
compounds 
405.101 C17H2209 M+CI]- Perilloside E Phenolic glycosides 11565 
157.0973 C7H12N202 M+H] (2Z)-1-Ethyl-2-(nitromethylidene) Pyrrolidines 1386 c 
pyrrolidine 
Up 590.1275 C28H27NO12 M+Na-2H]- Lactonamycin Anthracenes 18028 
regulation 504.0783 C17H23N5013 M-H]- UNPD127787 Carboxylic acids and 15377 
(22) derivatives 
552.1514 C28H27NO11 M-H]- Amidepsine B;FO 2942B Depsides and 16891 c¢ 
depsidones 
574.1329 C28H27NO11 M+Na-2H]- Amidepsine B Depsides and 17533, oc 
depsidones 
552.1513 C28H27NO11 M-H]- Amidepsine B;FO 2942B Depsides and 16889 c¢ 
depsidones 
243.066 C14H1204 M-H]- (2E,11Z)-Wyerone acid Fatty Acyls 5489 c 
313.2735 C19H3603 M+CH30H+H]+ = 18-oxo-nonadecanoic acid Fatty Acyls 4652 c 
441.0826 C22H18010 M-H]- (—)-Catechin 3-O-gallate Flavonoids 12960 c 
415.1454 C22H2408 M-H]- 8-Acetoxypinoresinol Furanoid lignans 11887 c¢ 
959.57 C51H84015 M-+Na]+ 1,2-Di-(9Z,12Z,15Z-octadecatrienoyl)-3- Glycerolipids 10139 c 
(Galactosyl-alpha-1-6-Galactosy]-beta- 
1)-glycerol 
452.2781 | C21H44NO7P M-H]- LysoPE(16:0/0:0) Glycerophospholipids 13357 
275.11 C12H1807 M+H]+ Xylomollin Lactones 3773 
540.3304 C28H43N704 M-H]- 1-[[(8S,9R)-6-[(2S)-1-hydroxy propan-2- Macrolactams 16593 


yl]-8-methyl-5-oxo-10-oxa-1,6,14,15- 
tetrazabicyclo[10.3.0]pentadeca-12,14- 
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Factors Down/Up m/z Formula Adduct ions Metabolite name Compounds ID QLs 
classification 
dien-9-yl]methyl]-1-methyl-3-[4-(1- 
piperidinyl)phenyl]urea 
699.3803 C38H56N2010 M-H]- Leucascandrolide A Macrolides and 20241 (a 
analogues 
340.2821 C18H39NO3 M+Na] Phytosphingosine Organonitrogen 5213 c 
compounds 
318.3002 C18H39NO3 M+H]+ Phytosphingosine Organonitrogen 4767 c 
compounds 
335.0948 €11H20010 M+Na] Serotinose Organooxygen 5079 c 
compounds 
215.0526 C7H1206 M+Na] Quinic acid Organooxygen 2540 b 
compounds 
487.215 C21H36011 M+Na] Linalool 3,6-oxide primeveroside Organooxygen 7756 b 
compounds 
353.2685 C21H3604 M+H]+ Sterebin K; (+)-Sterebin K Prenol lipids 5451 c 
531.0521 C15H22N2015P2 M-H]- [(2R,3R,4R)-3,4-dihydroxy-5-keto- Pyrimidine nucleotides 16264 ic 
tetrahydropyran-2-yl] [[(2R,3S,5R)-5- 
(2,4-diketo-5-methyl-pyrimidin-1-yl)-3- 
hydroxy-tetrahydrofuran-2-yl]methoxy- 
hydroxy-phosphoryl] hydrogen 
phosphate 
526.0599 C24H19N504S82 [M+Na-2H]- LSM-4864 Quinolines and 16120 c 


derivatives 


Note: Qualitative levels(‘‘a", Standards; “b" Laboratory self-built database-TMDB; “c", MS2). 


analysed (Fig. 5). Intriguing insights were obtained from the investiga- 
tion, which revealed impressive consistency in the changes of the 10 
enzymes under scrutiny (ANR, PPO, POD, FLS, PDAT, PDCT, LPCAT, 
pectase, ANS, and F3H) (Fig. 5b). The findings suggested that during the 
initial phase of fermentation, the production of enzymes was minimal, 
and a negligible number of enzymes were detected. However, from the 
second and third stages of fermentation, the amount of enzymes grad- 
ually increased. The highest enzyme activity was reached during the 
fourth and fifth stages of fermentation, when the activity of the enzymes 
exhibited an unprecedented surge. However, during the sixth stage of 
fermentation, the activity of the enzymes started to decrease noticeably. 
This intriguing pattern of enzyme activity closely aligns with the for- 
mation of premium Pu-erh tea. As per our understanding, the natural 
fermentation process is influenced by several environmental factors 
such as humidity, temperature, and the interplay of microorganisms, 
which collectively foster the production of enzymes, leading to a sig- 
nificant upsurge in the enzymes’ activity. 


3.4.2. Analysis and discussion of the biosynthesis of flavonoids and 
flavonols 

As displayed in Fig. 4a, the composition of Pu-erh tea was affected by 
the PF environment (microbial action and moist-heat action) mainly 
through two biosynthetic pathways: the phospholipid pathway and the 
flavone and flavonol biosynthetic pathway (Fig. 4a). As shown in Fig. 4a, 
upstream of the flavone, flavanols, and flavonoid biosynthesis biosyn- 
thetic pathway, the naringenin and dihydrokaempferol content signifi- 
cantly decreased during processing. However, we found that naringenin 
and dihydrokaempferol were significantly degraded by mainly micro- 
bial action; EGCG, EGC, and ECG were oxidatively degraded by moist 
heat (Fig. 5a). Downstream of the flavone, flavanols, and flavonoid 
biosynthesis biosynthetic pathway, we discovered synthesis of some 
taste-presenting substances (quercetin, kaempferol, GA, catechins 
oxidative polymers, and puerin I-IV). During long-term PF, the amounts 
of puerin I, puerin II, puerin II, puerin IV, quercetin, and kaempferol 
started to increase in the second stage of fermentation, reached their 
highest values in the fourth and fifth stages of fermentation, and then 
started to decrease (Fig. 5a). Combining this with the analysis of the 
relevant enzyme activities, we concluded that the activities of ANR, 
PPO, POD, FLS, and F3H were strongly positively correlated with the 
synthesis of metabolites such as GA, puerin I, puerin II, puerin III, puerin 
IV, quercetin, and kaempferol and negatively correlated with the 


synthesis of catechin monomers (EGC, EC, C, EGCG, ECG, and TC; 
Fig. 5c). The rate of synthesis of catechin monomers during the PF of Pu- 
erh tea was much lower than the rate of their oxidative degradation into 
puerin and GA. 

Analysis of volcano plots and Venn diagrams revealed that the me- 
tabolites quercetin, kaempferol, and GA were synthesised mainly 
through microbial effects (Li et al., 2022), whereas the metabolites 
dihydroisorhamnetin, apigenin 7-xyloside, apigenin 7-(3’-ace- 
tyl-6"-E-p-coumaroylglucoside), epicatechin quinone, and _ dihy- 
drokaempferol were mainly synthesised through moist-heat action. We 
also discovered that the metabolites GA, kaempferol, quercetin, and 
puerin I-IV and the catechin monomers (GC, C, EGCG, EC, and ECG) 
were synthesised through a combination of microbial action and 
moist-heat action (Table 1). Li et al. (2022) concluded that GA can be 
synthesised in large amounts under the action of humid heat (Fig. S1). 
This conclusion is consistent with the findings of the present study, but 
we found that the efficiency of GA synthesis under microbial action was 
significantly higher than that of GA synthesis under the effect of moist 
heat (Table 1). The oxidation products of catechins were rapidly 
degraded through oxidation under microbial action and further gener- 
ated substances such as theobromine and theaflavin. In the study con- 
ducted by Jiang et al. (2021), Jiang et al. (2022) and Jiang et al. (2023) 
substances such as N-ethyl-2-pyrrolidinone-substituted flavan-3-ols 
were found to promote the mellowing of tea infusions, whereas GA 
mainly promoted sweetness. Additionally, the review conducted by 
Wang et al. (2022) mentioned that GA is an aroma precursor for 
methoxyphenols, providing the basis for the stale aroma of Pu-erh tea. 
Amino acids and catechins are easily converted into puerarin under the 
action of moist heat. Amino acids, EGCG, and EGC produce 
flavan-alkaloids (N-ethyl-2-pyrrolidinone-substituted —flavan-3-ols). 
Upon analysing the quantitative levels of amino acids and 
catechin-like substances (as illustrated in Figures $1, $2 and S3), it was 
revealed that a significant portion of amino acids had undergone 
degradation as a result of microbial activity and exposure to moist heat 
during processing. Only a handful of amino acids, such as Gly, Orn, b-Ala 
and Cysthi, were found to be synthesised within the processing envi- 
ronment (Figs. S2 and S3). Studies have shown that the oxidation of 
theanine and catechin results in the production of puerin (flavan-alka- 
loids), which plays an important role in the quality of tea leaves (Zhang 
et al., 2013; Zhou, Zhang, Xu, & Yang, 2005). Therefore, we believe that 
this may be the reason for the bitterness and astringency of Pu-erh tea 
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Piling fermentation stage 3 


Trends in chemical composition 


Fig. 4. Key pathway analysis: (a) biosynthesis pathway of flavonoids and flavonol, and (b) possible lipid metabolite pathway during piling fermentation. Note: 
polyphenol oxidase (PPO), pectase, peroxidase (POD), anthocyanin reductase, anthocyanidin synthase (ANS), flavonol synthase (FLS), flavanone-3-hydroxylase 
(F3H), phosphatidylcholine diglyceride choline phosphotransferase (PDCT), lecithin acyltransferase (LPCAT), and diacylglycerol acyltransferase (PDAT). Different 
fermentation stages: Piling fermentation-1 (PF-1), Piling fermentation-2 (PF-2), Piling fermentation-3 (PF-3), Piling fermentation-4 (PF-4), Piling fermentation-5 (PF- 
5), Piling fermentation-6 (PF-6). The mini heat map from blue to red represents the activity from low to high. The green mini-bars represent the trend in the content 
of key chemical components during piling fermentation. (For interpretation of the references to colour in this figure legend, the reader is referred to the Web version 
of this article.) 


being reduced and its taste thickness and mellowness being increased. Li, 2015). The synthesis of phospholipid-like substances has been found 
to make a positive contribution to tea aroma (Li et al., 2023; Wang et al., 

3.4.3. Analysis and discussion of the lipid metabolite pathway 2022). To further investigate the mechanism of phospholipid-like sub- 
Lipids are an important component of tea leaves, accounting for stance changes, we tested related enzymes (PDAT, PDCT, LPCAT, and 


approximately 8% of the dry weight of tea leaves (Wan, 2003). Lipids pectase, as shown in Figs. 4b and 5b). Correlation analysis (Fig. 5d) 
play a pivotal role in the formation of tea quality and are considered one revealed that PDCT, LPCAT, pectase, and PDAT were positively corre- 
of the main pathways through which tea’s aroma forms (Ho, Zheng, & lated with the synthesis of lipids. Together with the above analysis (see 
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Fig. 5. Analysis of relevant metabolites and enzyme activities in key pathways under spontaneous fermentation. (a) Profile of key metabolite changes, (b) trends of 


different enzyme activities, (c) and (d) correlation analysis. 


Note: PF-CK, sundried green tea; SPPF, spontaneous piling fermentation. SPPF-processed samples [SPPF-1 (0 days), SPPF-2 (7 days), SPPF-3 (14 days), SPPF-4 (22 
days), SPPF-5 (32 days), SPPF-6 (38 days)]; MA, microbial action; MA-MHA, microbial action and moist-heat action; MHA, moist-heat action; polyphenol oxidase 
(PPO), pectase, peroxidase (POD), anthocyanin reductase, anthocyanidin synthase (ANS), flavonol synthase (FLS), flavanone-3-hydroxylase (F3H), phosphatidyl- 
choline diglyceride choline phosphotransferase (PDCT), lecithin acyltransferase (LPCAT), and diacylglycerol acyltransferase (PDAT). (For interpretation of the 
references to colour in this figure legend, the reader is referred to the Web version of this article.) 


3.3.3), we found that lysophospholipids were mainly regulated through 
the action of microorganisms. Li et al. found that during the fermenta- 
tion of Pu-erh tea, phospholipids were heavily degraded, whereas 
lysophospholipids first increased and then decreased in number. Li et al., 
2023 reported that fatty acid esters of hydroxy fatty acids could only be 
synthesised during the fermentation of Pu-erh tea, which is influenced 
by microbial effects. However, the study ignored the fact that microbial 
action, heat, and humidity during the fermentation of Pu-erh tea play 
crucial roles in the formation of the tea’s quality. Therefore, we believe 
that the results obtained by Li et al., 2023 have some limitations and 
cannot clearly indicate which factor causes the synthesis and degrada- 
tion of metabolites. In this study, our findings demonstrated that lyso- 
phospholipids could be synthesised only through the action of microbes 
and were unaffected by the action of moist heat. Therefore, we 
concluded that the synthesis of lysophospholipids is mainly a result of 
the microorganisms generated during the fermentation process that is, 
through the degradation and decomposition of tea leaves. Secondly, we 
will focus more on the origin of microorganisms and the dynamic evo- 
lution process of microorganisms in the processing environment of 
Pu-erh tea in our future research, and thus investigate the mechanism of 
tea quality formation. 


4. Conclusion 


In this study, raw SDGT was used for the first time to investigate the 
effects of microbial action and moist-heat action on metabolite regula- 
tion during the PF of Pu-erh tea. The sensory evaluation revealed that 
the colour of a Pu-erh tea infusion is brownish-red and that the infusion 
has a stale, fungal stale, and pure aroma as well as a stale and mellow 
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taste when the leaves are left to ferment naturally. When tea leaves 
undergo STPF, the infusion is orange-red and dull-red in colour, has a 
woody and sweet aroma, and has a sweet and mellow taste. We screened 
108 metabolites by performing histological analysis. A total of 25 me- 
tabolites were significantly up-regulated under microbial action, 
whereas 8 were significantly downregulated; 22 metabolites were 
upregulated and 17 were downregulated by moist-heat action; and 
finally, 10 metabolites were downregulated and 26 were upregulated by 
the combined actions of microbes and moist heat. Combining enzyme 
activity and pathway analysis, we determined that the enzymes pro- 
duced by microorganisms play a crucial role in expediting oxidative 
degradation through the flavonoid metabolic pathway and lipid 
metabolite pathway, ultimately leading to development of the distinct 
quality of Pu-erh tea. This study clearly discovered that the lipid 
metabolite pathway is regulated by the action of microorganisms, thus 
influencing the aroma formation of Pu-erh tea. In conclusion, we found 
that these two pathways are crucial in the formation of the stale aroma 
and stale-mellow taste of Pu-erh tea. This study complements our pre- 
vious research and provides a theoretical basis for identification of the 
quality formation of Pu-erh tea. 
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